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Open 



Microbial fuel cells (MFCs) are often inoculated from a single wastewater source. The extent that the 
inoculum affects community development or power production is unknown. The stable anodic 
microbial communities in MFCs were examined using three inocula: a wastewater treatment plant 
sample known to produce consistent power densities, a second wastewater treatment plant sample, 
and an anaerobic bog sediment. The bog-inoculated MFCs initially produced higher power densities 
than the wastewater-inoculated MFCs, but after 20 cycles all MFCs on average converged to similar 
voltages (470±20mV) and maximum power densities (590±170mWm 2 ). The power output from 
replicate bog-inoculated MFCs was not significantly different, but one wastewater-inoculated MFC 
(UAJA3 (UAJA, University Area Joint Authority Wastewater Treatment Plant)) produced substantially 
less power. Denaturing gradient gel electrophoresis profiling showed a stable exoelectrogenic 
biofilm community in all samples after 11 cycles. After 16 cycles the predominance of Geobacter 
spp. in anode communities was identified using 16S rRNA gene clone libraries (58 ±10%), 
fluorescent in-situ hybridization (FISH) (63 ±6%) and pyrosequencing (81 ±4%). While the clone 
library analysis for the underperforming UAJA3 had a significantly lower percentage of Geobacter 
spp. sequences (36%), suggesting that a predominance of this microbe was needed for convergent 
power densities, the lower percentage of this species was not verified by FISH or pyrosequencing 
analyses. These results show that the predominance of Geobacter spp. in acetate-fed systems was 
consistent with good MFC performance and independent of the inoculum source. 
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Introduction 

In a microbial fuel cell (MFC), exoelectrogenic 
bacteria on the anode degrade organic matter and 
produce electrical current. Exoelectrogens occur 
widely in nature (Logan, 2009), and MFCs have 
been successfully inoculated using pure cultures 
(Nevin et ah, 2008; Ishii et ah, 2008b; Watson and 
Logan, 2010), aerobic and anaerobic wastewater 
streams (Min et ah, 2005; Ki et ah, 2008), as well 
as environmental samples such as seawater, rivers 
and soil (Bond et ah, 2002; Ieropoulos et ah, 2010; 
Jiang et ah, 2010). 

The role of different inoculum sources on the 
communities that develop and the associated MFC 
performances under otherwise identical experi- 
mental conditions has not been well studied. 
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Holmes et ah (2004) examined anodic communities 
and compared reactor performance of sediment-type 
MFCs inoculated with samples from salt-water 
marshes, freshwater and marine sediments. They 
found that all anodic communities were dominated 
by different orders of Deltaproteobacteria, indicating 
that this class of bacteria was linked to exoelectro- 
genic activity. MFCs inoculated with wastewater 
have produced more power than a natural environ- 
ment inoculum in some cases, but not in others. 
Anaerobic and aerobic wastewater effluent inocula 
were compared with a river water inoculum in a 
continuous-flow, acetate-fed MFC (Ieropoulos et ah, 
2010). The aerobic effluent (1.16 mWm" 2 ) and river 
water (1.14 mWm" 2 ) inocula produced a higher 
maximum power than the anaerobic effluent inocu- 
lum (0.89 mWm" 2 ). In contrast, a wastewater 
inoculum produced a higher open circuit 
voltage (785 mV) than a soil inoculum (684 mV) 
(Jiang et ah, 2010) when fed acetate, but power 
densities were not reported. The communities 
that developed in these wastewater studies were 
not reported. 
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Another issue that arises in many MFC studies is a 
lack of replicate reactors. It has been noted that 
while replication is important to achieve statistical 
confidence in molecular community analyses, 
~60% of microbial diversity studies do not include 
replicate analyses (Prosser, 2010). Several MFC 
studies have been conducted using only a single 
reactor (Bond et ah, 2002; Aelterman et ah, 2006; 
Jong et ah, 2006; Reimers et ah, 2007; Xing et ah, 
2008b). The lack of replicates makes it difficult to 
determine if the reported microbial communities 
and reactor performances would be reproducible. 
Fortunately, there seems to be more of a trend in 
recent studies to use replicates, with most studies 
reporting on duplicate (Holmes et ah, 2004; 
Rismani-Yazdi et ah, 2011; Sasaki et ah, 2011) or 
triplicate reactors (Faimali et ah, 2010; Nelson et al., 
2010; Jung and Regan, 2011). It is particularly 
important to establish better linkages between 
reactor performance and microbial communities 
because other factors affect reactor performance, 
such as differing reactor internal resistances and 
architectures. 

Comparisons of microbial communities are also 
complicated by the fact that many studies use only 
one molecular technique to analyze the anodic 
community. The two most commonly used commu- 
nity analysis techniques in MFC studies have been 
denaturing gradient gel electrophoresis (DGGE) 
(Rismani-Yazdi et al., 2007; Erable et al., 2009; 
Zang et al., 2010; Jung and Regan, 2011) and clone 
libraries (Ishii et al., 2008a; Ki et al., 2008; Kiely 
et al, 2011a; Kiely et al., 2011b; Sun et al., 2010) 
using PCR amplicons of 16S rRNA genes. Pyrose- 
quencing can be used to increase the sensitivity of 
community analysis compared with DGGE or clone 
libraries, allowing identification of less-abundant 
(but possibly still important) members of the com- 
munity, while eliminating cloning bias (Lee et al., 
2010). However, all of these post-PCR methods are 
potentially subject to quantitative bias associated 
with differential amplification efficiency and rRNA 
operon copy number. The use of multiple analysis 
techniques may provide greater confidence on the 
composition of the community. For example, fluor- 
escent in-situ hybridization (FISH) is a useful non- 
PCR based technique that can be used to corroborate 
the findings of other molecular techniques (Kiely 
et al., 2011c), but it has rarely been used in MFC 
studies. 

In this study we used four different molecular 
analysis techniques [16S rRNA gene-targeted DGGE, 
clone libraries and pyrosequencing and FISH) to 
study the anodic communities that develop in 
MFCs. Three different inocula were studied in 
triplicate reactors over time. Through these ana- 
lyses, we were able to link reactor performance to 
shifts in community composition and predomi- 
nance of important exoelectrogenic species in MFCs 
that were otherwise identical in architecture and 
chemical solutions. 



Materials and methods 

Reactor configuration, inoculation and electrochemical 
monitoring 

Single-chamber air-cathode MFCs (28 ml) were con- 
structed from Lexan cubes to form a 4-cm long 
cylindrical chamber 3 cm in diameter (Cheng and 
Logan, 2007). Anodes were made from a graphite 
fiber brush (2.5 cm diameter, 2.5 cm long) (Logan 
et ah, 2007) (PANEX33 160 K, ZOLTEK, St Louis, 
MO, USA) and cathodes were carbon cloth (type 
B-1B, E-TEK) containing 0.5 mg cm" 2 Pt (Cheng 
et ah, 2006). Reactors were kept in a constant 30 °C 
chamber and fed a medium consisting of 1 g 1 ~ 1 
sodium acetate and Wolfe's vitamins (5 mil -1 ) and 
minerals (12.5 ml!- 1 ) (Balch et ah, 1979) in 50 mM 
phosphate buffer (Na 2 HP0 4 , 4.56 gl" 1 ; NaH 2 P0 4 
monohydrate, 2.45 gl" 1 ; NH 4 C1, 0.31 gl" 1 ; KC1, 
0.13 gl" 1 ) (Xing et ah, 2009; Kiely et ah, 2011b). 
There was no attempt made to prevent light 
intrusion into the reactors. 

Reactors (triplicate) were inoculated using samples 
from three different sources: the PSU (Pennsylvania 
State University) Wastewater Treatment Plant primary 
clarifier effluent, which has been used in a number 
of MFC studies with relatively good and consistent 
performance (Liu et ah, 2005; Kiely et ah, 2011b); 
University Area Joint Authority Wastewater Treat- 
ment Plant (UAJA) primary clarifier effluent; 
and sediment from a stagnant bog in the Black 
Moshannon Lake State Park, PA, USA (Bog). The 
UAJA plant is much larger, and serves more rural 
communities, compared with the PSU plant that 
treats 'city' wastewater. Wastewater reactors were 
inoculated in a 1:1 ratio of wastewater to medium. 
The Bog reactors were inoculated using a 1:3 
ratio of sediment to medium to maintain similar 
viscosity to the wastewater inocula. All reactors 
were then fed every 3 days (with 3 days being 
equivalent to one cycle), with complete medium 
replacement. 

Voltage [E) was measured every 20min across an 
external resistor [R = 1 kQ) using a data collection 
system (Keithley, Cleveland, OH, USA). Power 
[P = E 2 /R) and current [I=E/R] were normalized by 
the cathode projected surface area (7.1cm 2 ). Polar- 
ization curves were generated using the single-cycle 
method (Watson and Logan, 2011) by varying 
external resistances every 20min from 10 kQ to 
50 Q. All statistical calculations using a £-test were 
based on a = 0.05. The Coulombic efficiency in 
similar fixed-resistance systems has been well 
studied (~20% for an external resistance of 
1000 Q) and was not measured in these experiments 
(Zhang et ah, 2009). 



DNA extraction 

Anode samples were taken during cycles 1, 3, 4, 5, 
11 and 16 by cutting ~0.2 g (wet weight; corre- 
sponding to ~2% of the anode) of the anode fibers 
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using flame-sterilized scissors. The samples 
were taken mid-cycle during peak electrochemical 
activity. Anode removal would not affect perfor- 
mance as recent tests with identical reactors have 
shown that up to 75% of the anode brush can be 
removed without adversely affecting voltage or 
power generation (Hutchinson et al, 2011). DNA 
was extracted using the PowerSoil DNA Isolation 
Kit (MO-BIO Laboratories, Carlsbad, CA, USA) 
according to the manufacturer's instructions. DNA 
samples from all cycles were used for cluster 
analysis of DGGE profiles via Principal Component 
Analysis to screen for convergence of the anodic 
communities from the different inocula. The con- 
vergent anode samples obtained during cycle 16 
were then further analyzed by FISH and used for the 
clone library and pyrosequencing analyses. DNA 
samples were stored at - 20 °C before analysis. 

1 6S rRNA gene clone libraries 

16S rRNA gene clone libraries were generated for 
the anode communities present during cycle 16 as 
previously described (Kiely et al, 2011b), but with 
different primers. Extracted DNA was amplified via 
PCR (Zuo et al, 2008) with a static annealing 
temperature of 55 °C to obtain fragment lengths of 
960 bp using bacterial primers 530F (5'-GTCCCAGC 
MGCCGCGG-3 / ) and 1490R (5'-GGTTACCTTGTTAC 
GACTT-3') (Wani et al, 2006). NCBI BLAST (blastn 
algorithm) was then used to analyze the DNA 
sequences using the NCBI Genomes database 
(Altschul et al, 1990). 

Shannon diversity indices, Good's coverage 
values and rarefaction curves were calculated using 
the Mothur computer software (Schloss et ah, 2009). 
Sequences were first aligned in MEGA 4 (Tamura 
et ah, 2007). Rarefaction curves are an average 
composed of 1000 randomizations. Approximately 
90 clones were analyzed for each sample and 
deposited into GenBank under accession numbers 
JF817397-JF818119. 



Pyrosequencing and analysis 

Nine additional bacterial phylogenetic libraries 
were produced from the cycle 16 anode samples 
by pyrosequencing utilizing the GS FLX sequencer 
and Titanium series chemistry (454 Life Sciences, 
Branford, CT, USA). Universal bacterial PCR primers 
343F S'-TACGGRAGGCAGCAG-S' and 926R 5'-CCG 
TCAATTYYTTTRAGTTT-3' (Liu et al, 2007; Wang 
and Qian, 2009) were used and also included 
pyrosequencing adapters, keys and multiplex iden- 
tifiers. PCR was conducted in 50 |il volumes contain- 
ing 1 x PCR master mix (Roche Applied Science, 
Indianapolis, IN, USA), 0.4 jim of each primer and 
2 |il of DNA template (5 ng |il " 1 ). PCR was performed 
at the following cycling conditions: initial denatura- 
tion at 94 °C for 5min, and 25 cycles of 94 °C 
denaturation for 15 s, annealing at 56 °C for 45 s 



and a 1-min extension at 72 °C, followed by a final 
extension at 72 °C for 8 min. Four reactions were 
conducted for each sample and amplicons were 
combined before removing salts and unincorporated 
primers using a Qiagen MinElute PCR purification 
kit (Qiagen Inc., Valencia, CA, USA) (Bibby et al, 
2010). Amplicons were visualized on a 1.2% agarose 
gel and extracted if smearing was observed using a 
Qiagen MinElute gel extraction kit (Qiagen Inc.). 
Each purified DNA amplicon was pooled at equal 
concentrations for sequencing. Before sequencing, 
DNA concentration and purity were confirmed on 
a 2100 Bioanalyzer (Agilent Technologies, Santa 
Clara, CA, USA). Sequencing was performed using 
an emPCR Lib-A Kit (Roche Applied Science) at 
the Yale Center for Genome Analysis. Reads were 
subjected to quality control at the machine based 
on run-produced Phred scores. Keypass, dots and 
mixed filters assessed the quality of the whole read 
while the quality of read ends was checked by signal 
intensity and primer filters. 

Quantitative sequence analysis was performed 
using programs in the Quantitative Insights Into 
Microbial Ecology (QIIME) tool box (Caporaso et al., 
2010). Within QIIME, sequences were sorted, and 
primer and multiplex identifier sequences were 
trimmed. Denoising was conducted using Titanium 
Pyronoise software (Quince et ah, 2009; Caporaso 
et ah, 2010; Quince et ah, 2011). Similar sequences 
were clustered into operational taxonomic units 
based on 97% identity, aligned with the Greengenes 
core set (DeSantis et al., 2006), and taxonomy 
was assigned using the RDP (ribosomal database 
project) classifier (Cole et al, 2009). Within QIIME, 
rarefaction curves using the denoised data were 
generated using the alpha_rarefaction.py workflow 
script with OTUs clustered at 97% identity. The 
unprocessed DNA sequences obtained in this study 
have been deposited in the MG-RAST (Meyer et al., 
2009) archive under accession number 4465374.3. 
Multiplex identifier information is included in 
Supplementary Table Si. 



DGGE analysis 

Fragments of 16S rRNA genes from cycles 1 and 16 
were amplified via PCR using universal primers 
968F (5 / -AACGCGAAGAACCTTAC-3 / ) and 1401R 
(5 / -CGGTGTGTACAAGACCC-3 / ). A GC clamp (5'-CG 
CCCGCCGCGCCCCGCGCCCGGCCCGCCGCCCCCGC 
CCC-3') was attached to 968F (Ren et al, 2007). 
Samples were amplified with a G-Storm 1 Thermo- 
cycler (GRI Laboratories, Braintree, Essex, UK) using 
a previously described temperature profile (Xing 
et al, 2010) and stored at 4 °C. DGGE gels were 
prepared as previously described (Xing et al, 
2008a). Electrophoresis was conducted for 12 h at 
60 °C and 70 V in a 0.5 x TAE buffer solution to 
achieve effective band separation. Gels were then 
silver stained as previously described (Bassam 
et al, 1991). 
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The most visually prominent bands from each gel 
were excised using sterilized pipette tips, crushed in 
50 |il of TE buffer (10 mM Tris-HCl, ImM EDTA (pH 
8.0)) and eluted overnight at 4 °C. The eluted DNA 
was used as the template for PCR as described above 
excluding the GC clamp. The PCR products were 
purified using a PCR cleanup kit (Qiagen) followed 
by ligation of the DNA fragments into a TOPO 2.1 
vector (Invitrogen, Carlsbad, CA, USA) and trans- 
formation into chemically competent E. coli cells. 
Plasmids were extracted (96-well Plasmid Extrac- 
tion Kit, Omega Bio-Tek, Norcross, GA, USA) from 
eight random colonies per band and the insert DNA 
was sequenced (ABI 3730XL, Applied Biosystems, 
Carlsbad, CA, USA) to phylogenetically identify the 
bands and inspect for comigration of bands derived 
from different templates. 

PCA of the DGGE profiles was conducted using 
two methods in order to determine the convergence 
of the anode communities. Both methods used 
presence/absence of bands and scoring band inten- 
sities (scale from 0-5). One method used a computer 
software program (GelCompar II, Applied-Maths, 
Austin, TX, USA), while the other used visual 
observation (Fry et al., 2006). PCA plots were 
generated using Minitab 16 (Minitab, State College, 
PA, USA) to determine the similarity of the 
replicates at different times and to evaluate temporal 
community shifts in the reactors. 



FISH 

The final anode sample taken during cycle 16 was 
analyzed using FISH for comparison with clone 
library and pyrosequencing results. Samples were 
immediately fixed using 4% paraformaldehyde 
for 3 h at 4 °C, washed in sterile phosphate buffer 
(FISH-PBS; lOmM NaH 2 P0 4 , 130mM NaCl, pH 7.2-7.4), 
and stored at -20 °C in 1:1 PBS to ethanol solution. 
A G. sulfurreducens specific probe 5 r -labeled with 
Alexa Fluor 594 (GE02; 5'-GAAGACAGGAGGCCCG 
AAA-3') (Invitrogen) along with 2 unlabeled helper 
probes (HGE02-2; 5'-CTAATGGTACGCGGACTCA 
TCC-3' and HGE02-1; S'-GTCCCCCCCTTTTCCCGC 
AAG-3') (Richter et al, 2007) were each added to a 
final concentration of 10ng|il -1 and hybridized in 
1m NaCl solution containing 0% formamide for 
1.5 h at 46 °C (Kiely et al, 2010). Samples were then 
washed in a 1-m NaCl solution at 48 °C and stained 
with Slow Fade Gold containing DAPI (Invitrogen) 
to identify all bacterial and archaeal cells. Visualiza- 
tion of DAPI and probe-conferred fluorescence was 
then performed on an Olympus BX61 epifluorescent 
microscope. The percentage of G. sulfurreducens in 
each reactor sample was determined by counting at 
least 1000 total cells from 10 different fields. No 
probe and nonsense probe controls were included 
to examine for autofluorescence and non-specific 
binding of the dye, respectively. Verification of the 
FISH probe specificity was performed using probe- 
Check (Loy et al., 2008) and through negative 



hybridizations with two closely related bacteria 
previously found in acetate-fed MFCs [Pelobacter 
propionicus, G. metallireducens) (data not shown). 
Each species tested had six mismatches from the 
probe used in this study. 

Results 

Electrochemical performance 

The bog-inoculated reactors initially produced 
higher maximum voltages on average than the 
wastewater-inoculated MFCs, but after 20 cycles, 
all reactors converged to a maximum voltage output 
of 470±20mV (Figure 1). Based on polarization 
tests, maximum power densities on cycle 16 were 
not significantly different for eight of the nine 
reactors (P>0.05), but the UAJA3 replicate pro- 
duced significantly (P= 0.005, Mest) less voltage 
(data not shown) and power (354 mWm" 2 ; Figure 2) 
than the other samples. The maximum power 
density was 590±110mWm" 2 based on all nine 
reactors, and 620 ± 60 mWm" 2 excluding the UAJA3 
replicate (n = 8). 

16S rRNA gene clone library and pyrosequencing 
analysis 

Based on the 16S rRNA gene clone libraries, during 
cycle 16 all reactors were found to be predominated 
by clones most similar to Geobacter spp. (>95% 
identity to Geobacter spp. with most sequences 
having >97% identity to G. sulfurreducens). On 
average, 61 + 6% (27 = 8, excluding UAJA3) of the 
clones were similar to Geobacter spp., except for 
replicate UAJA3, which had 36%. Clones distantly 
related to Bacteroidetes (89% identity or lower) 
were the second most abundant taxonomic group in 
all of the anode communities (11 + 6%, n = 8), but 
were 36% of the clone library for the community in 
the UAJA3 replicate. Other sequences identified in 
the clone analysis were <95% similar to members 
of Betaproteobacteria (5 + 3%, n = 9) and Alphapro- 
teobacteria (2 + 2%, n = 9) (Figure 3). 



0.6 




Cycle Number 

Figure 1 Average maximum voltage [n = 3) in each cycle for each 
inoculum source (B=bog, P = PSU and U = UAJA). 
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Figure 2 Power density curves generated from cycle 16 showing 
that UAJA replicate 3 produced far less power than all other 
reactors, including other UAJA replicates. Legend order mimics 
the order of the curves based on max power density. 



A higher relative abundance of sequences with 
similarity to Deltaproteobacteria (81 + 4%) (Figure 4) 
was identified by pyrosequencing compared with 16S 
rRNA gene clone libraries. Bacteroidetes (10 + 5%) 
and Firmicutes (8 + 3%) were the other two most 
prominent phyla present. The Deltaproteobacteria 
sequences obtained by pyrosequencing were domi- 
nated by the closely related genera Geobacter and 
Trichlorobacter, as identified by the RDP II classifier 
(Cole et ah, 2005). It has been proposed that 
Trichlorobacter be reclassified as Geobacter (Nevin 
et ah, 2007), so the two genera are grouped together in 
this analysis. The main difference between the two 
results is that pyrosequencing analysis did not 
confirm the clone library result that showed a lower 
abundance of clones with sequence similarity to G. 
sulfurreducens in the UAJA3 replicate. 
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Figure 3 Abundance of clones with sequence identity to dominant classes in each reactor from cycle 16. UAJA3 had the lowest number 
of clones showing similarity to Deltaproteobacteria. The group labeled as other consists of sequences with <95% identity. Insert: the 
majority of the Deltaproteobacteria were Geobacter sulfurreducens. 
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Figure 4 Community analysis results obtained from pyrosequencing. No significant difference in the amount of Proteobacteria was seen 
in any reactor. Insert: the Deltaproteobacteria was made up almost entirely of Geobacter and Trichlorobacter (95%). Other minor genera 
were detected by pyrosequencing that were not detected in the clone libraries. 
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Table 1 Shannon Diversity Indices and Good's coverage values 
calculated from 1 6S rRNA gene clone library and pyrosequencing 
data taken during cycle 16 



16S rRNA clone libraries 



Pyrosequencing 





No. of 


Shannon 


Good's 


No. of 


Shannon 




clnrtp^ 


d iversi ty 




O C7 \J Li C71 O 


iui v C7i ji ty 


Bog 1 


93 


1.57 


86 


3413 


3.82 


Bog 2 


92 


1.94 


80 


1757 


4.04 


Bog 3 


92 


1.45 


89 


1698 


4.07 


Bog- 




1.65 


85 




3.98 


average 












PSU 1 


86 


1.86 


81 


3073 


5.07 


PSU 2 


94 


1.53 


84 


2675 


4.36 


PSU 3 


95 


1.99 


82 


2745 


4.32 


PSU- 




1.79 


82 




4.58 


average 












UAJA 1 


95 


1.90 


79 


1022 


4.39 


UAJA 2 


88 


1.72 


79 


1043 


4.25 


UAJA 3 


91 


2.07 


81 


3595 


4.69 


UAJA- 




1.90 


80 




4.44 


average 
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Figure 5 PCA plot of the dominant (darkest) bands from the 
DGGE gels representing the communities at cycle 1 and cycle 16. 
This plot contrasts the initial community of wastewater-inocu- 
lated reactors (P = PSU and U = UAJA) with the bog-inoculated 
reactors (B). Over the duration of the experiment all communities 
converged irrespective of inoculum source. 



Abbreviations: PSU, Pennsylvania State University; UAJA, University 
Area Joint Authority Wastewater Treatment Plant. 



There was no difference in diversity among the 
anode communities derived from the different 
inocula at cycle 16 based on the Shannon Diversity 
Index of each reactor generated from clone library 
data (1.78 ±0.22, n = 9, P>0.05). Good's Coverage 
values indicated that on average 82% of the total 
species richness was accounted for in each commu- 
nity (Pommier et ah, 2007) (Table 1). Rarefaction 
curves generated using both Mothur and QIIME 
(454 sequences) also suggest that there is little 
difference between the diversity of each of the 
reactors. The rarefaction curves using both data sets 
approached, but did not achieve, a plateau despite 
using more than 1500 sequences for the pyrosequen- 
cing data set compared with ~100 sequences using 
the 16S rRNA gene clone libraries. While this might 
suggest that the diversity of the bacterial commu- 
nities was not fully captured here, it is clear that 
the predominance of Geobacter was shown using 
both approaches. 



Community shifts with DGGE 

Analysis of predominant bands in DGGE gels using 
PCA showed convergence of the final microbial 
communities at cycle 16 (Supplementary Figure Si). 
During the first cycle, the Bog results clustered 
separately (quadrant 1) from the two wastewater 
inoculum sources (quadrant 4) (Figure 5). By the end 
of the study, all of the data clustered in quadrant 2. 
The computer-generated PCA (Warmink and van 
Elsas, 2008) produced the same results as manual 
scoring, indicating no bias in the manual scoring 
method (Supplementary Figure S2). 



2 1 




Figure 6 DDGE gel from cycle 16 showing bands that were 
excised and identified. Gels were used in the PCA to monitor 
shifts in overall microbial community throughout the experiment. 
(B=bog, P = PSU and U = UAJA). Asterisks denote bands that 
exhibited comigration during analysis. 



Seven of the visually darkest bands in the 
DGGE gels during cycle 16 were excised and 
sequenced, assuming that these bands would 
represent the predominant species in the com- 
munities (Figure 6). None of the dark bands 
analyzed (bands 2, 3 and 4) showed sequence 
identity to Geobacter spp., which was shown to be 
dominant in the 16S rRNA gene clone library and 
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Table 2 DGGE band identification 





Closest known isolate 


Identity 


Found in 






(%) 


16S rRNA gene 








clone library 


Band 1 


Geobacter sulfurreducens 


99 


Yes 


Band 2 


Alkaliphilus oremlandii 


89 


Yes 


Band 3 


Pseudomonas fluorescens 


99 


Yes 


Band 4 


Clostridium sticklandii 


89 


Yes 


Band 5 


Candidatus Protochlamydia 


89 


Yes 




amoebophila 






Band 6 


Rhodothermus marinus 


83 


No 




Bartonella quintana 


81 


Yes 


Band 7 


Exiguobacterium sp. 


88 


No 




Clostridium tetani 


84 


Yes 



Abbreviation: DGGE, denaturing gradient gel electrophoresis 




Figure 7 Percentage of G. suljurreducens in each reactor 
determined by FISH analysis, based on 10 images taken from 
different locations of each anode sample. 



pyrosequencing analysis. Analysis of additional 
bands (1 and 5) revealed the presence of bacteria 
that were most closely related to the known isolates 
Geobacter spp. (95% identity) and 'Candidatus 
Protochlamydia sp.' (89%). Other bands sequenced 
were most closely related to the known Firmicutes 
isolates Clostridium sp. (89%) and Alkaliphilus sp. 
(89%), and the Gammaproteobacterium Pseudo- 
monas sp. (99%) (Table 2). Cloning and sequencing 
of these prominent DGGE bands did not show 
evidence of DNA comigration. Other sequenced 
bands (6 and 7) did exhibit comigration, but did 
not include Geobacter. 



FISH 

FISH confirmed that a large portion of the 
biofilm communities was G. sulfurreducens. On 
the basis of probe-positive cell counts, 63 ± 6% of 
DAPI-stained cells were G. sulfurreducens [n = 30 
per sample) (Figure 7). FISH analysis did not 
confirm the clone library results that the UAJA3 
sample had a significantly smaller population of 
G. sulfurreducens. 



Discussion 

Voltage and power comparisons 

After a long acclimation time (20 cycles), the 
voltages of all the reactors converged to a similar 
voltage (470±20mV) and power density (590 + 
HOmWm" 2 ) regardless of inoculum source. This 
conclusion would not be different based on any 
combination of duplicates. The average voltages and 
maximum power calculated using the two most 
disparate replicates from the UAJA-inoculated reac- 
tors produced 450±34mV and 490 ± 190 mWm 2 , 
compared with 453±25mV and 510±140mWm 2 
based on triplicate reactors. These values for voltage 
(P=0.90) and power (P=0.92) are not significantly 
different when based on duplicates or triplicates. 
Power production in MFCs with flat or brush anodes 
is limited by the cathode (Fan et ah, 2008; Rismani- 
Yazdi et ah, 2008). MFCs with brush anodes have 
been shown to produce lower internal resistances 
and larger power densities than reactors with flat 
anode materials (Logan et ah, 2007), and therefore 
the systems used here were all cathode-limited in 
power generation. Even under the conditions where 
the anode did not limit power generation, the 
average performance of the reactors was indepen- 
dent of the initial inoculum once all reactors 
reached steady conditions. Given that UAJA3 
produced significantly less power than the other 
replicates by the end of the study, our conclusions 
could have been different based on the analysis of a 
single reactor. Although the Bog sample initially 
produced a significantly higher voltage than the 
wastewater-inoculated reactors, and these MFCs 
reached their maximum voltage after only three 
cycles, they eventually became similar in perfor- 
mance to the other inocula after 17 more cycles. This 
more rapid acclimation was likely due to consis- 
tently anaerobic conditions for microbial commu- 
nities in the bog compared with a wastewater 
effluent. 

Wastewater from the PSU treatment plant has 
previously been used in a large number of studies 
(Liu et ah, 2005; Cheng and Logan, 2007; Logan 
et ah, 2007; Xing et ah, 2010; Kiely et ah, 2011b) 
with results generally consistent with those reported 
here. Maximum voltages produced here with this 
and the other inocula (470±20mV) were generally 
consistent with those in several previous studies 
(500±40mV, 77 = 5) where 50 mM PBS, lgl" 1 
sodium acetate and single-chamber air-cathode 
MFCs were also used (Liu et ah, 2005; Cheng and 
Logan, 2007; Logan et ah, 2007; Xing et ah, 2010; 
Kiely et ah, 2011b). Maximum power densities have 
been slightly larger in some cases (704 ± 30 mW m " 2 , 
after the cathodic biofilm was removed) (Kiely et ah, 
2011b) and smaller in others (506 mWm" 2 ) (Liu 
et ah, 2005) than obtained here (590 ± llOmWrn" 2 ) 
with the same buffer. In some studies, power 
densities were only reported at higher buffer 
concentrations, precluding direct comparisons 
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(Cheng and Logan, 2007; Logan et al, 2007; Xing 
et al., 2010). Differences in maximum power for 
otherwise identical conditions could be due to 
methods used to obtain polarization data, which 
can result in higher power densities owing to Type 
M power overshoot (Watson and Logan, 2011), the 
use of different resistors or acclimation times before 
recording voltages; biofilms on the cathode, which 
can reduce performance; or better intrinsic perfor- 
mance of the cathodes leading to overall higher 
power densities. Some power density curves exam- 
ined here also exhibited Type D power overshoot, 
where the power density curve doubles back instead 
of reaching higher current densities. We now know 
that this type of overshoot results in an under- 
estimation of maximum power densities (Hong 
et al, 2011; Watson and Logan, 2011). Type D 
overshoot arises from the inability of the bacterial 
community to transfer electrons fast enough to the 
anode when the external resistance is lowered 
(Watson and Logan, 2011). Power overshoot implies 
that the capacity for these reactors to function at 
higher current densities would have required addi- 
tional acclimation to low resistances, which was not 
done here (Watson and Logan, 2011). 



Comparison of the community composition using 
different molecular techniques 

Analysis of communities from the different inocula 
showed similar amounts of G. sulfurreducens based 
on 16S rRNA gene clone library analysis (61 + 6%) 
and FISH (63 ± 6%), but a higher amount based on 
pyrosequencing data (81 + 4%). A similar conclu- 
sion of the predominance of Geobacter in acetate-fed 
MFC anode communities has been previously 
reported (Chae et al, 2009; Kan et al, 2011). In this 
study, acetate positively selected for Geobacter, 
which had a key role in its dominance of the 
communities in this study. Although Geobacter is 
seen to dominate in acetate-fed reactors, it is also 
able to grow using lactate (Call and Logan, 2011) and 
syntrophically using formic acid (Sun et al., 2012) 
showing that Geobacter is prevalent in systems 
using several different substrates. PC A of DGGE 
band patterns indicated a convergence of the 
communities to a similar composition after 16 
cycles, but sequencing of the darkest bands did not 
reveal a predominance of G. sulfurreducens or any 
other specific genus in the community. One possible 
cause for the differences among some of these 
techniques is that different primer sequences were 
used for each technique because of methodological 
constraints such as amplicon length for full dena- 
turation to give effective separation via DGGE (Wu 
et al., 1998) and sequencing length for pyrosequen- 
cing. The primer sets used here were chosen based 
on their common use in previous studies (Table 3). 
The possible effect these primers had on the 
resulting communities was analyzed with the RDP 
probe match function for sequences > 1200 bp with 
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no mismatches (Cole et al., 2009). In this analysis, 
the primers used here were compared with other 
bacterial primers reported in the literature for their 
identity to Deltaproteobacteria, Geobacter and Bac- 
teroidia, three of the most prevalent taxa in the clone 
library and pyrosequencing analyses (Table 3). It 
should be noted that low values for reported primers 
that target the extreme ends of the 16S rRNA gene 
could result from relatively fewer sequences cover- 
ing these regions in the RDP database. There also 
could be an inflation of matches for a particular 
primer if some sequences in the database were 
originally derived from the use of that primer and 
the primer sequence was not trimmed from the 
deposited sequence. The match scores for the DGGE 
primers used here are lower than those obtained for 
other common DGGE, clone library and pyrosequen- 
cing primers for all three taxa. This could in part 
explain the absence of Geobacter as a predominant 
community member in the DGGE analyses. Addi- 
tionally, the greater specificity of the pyrosequen- 
cing primers to Geobacter compared with the clone 
libraries primers could partially explain the differ- 
ence in Geobacter abundance. However, both primer 
sets have differences in targeting Bacteroidia, yet 
clone libraries (11 + 6%) and pyrosequencing 
(10 + 5%) gave similar results for sequences in this 
class. Despite these differences, in the final analysis, 
relatively similar responses were obtained for the 
different wastewater inocula based on both commu- 
nity analysis and time to reach maximum power 
densities. This indicates that these wastewater 
sources were equivalent as inocula. However, they 
were less effective in terms of achieving rapid 
acclimation of the reactors compared with the bog 
inoculum. 

Pyrosequencing analysis provided more depth of 
coverage to the community analysis than the 16S 
clone library analysis, which allowed a more 
detailed investigation of the bacteria present at low 
abundances in the anode community, indicating the 
presence of Acidobacteria, which includes the 
known exoelectrogen, Geothrix fermentans (Bond 
and Lovley, 2005) and Actinobacteria phyla, which 
were not identified using 16S rRNA gene clone 
libraries. The role that these members of the under- 
represented phyla could have in electricity genera- 
tion, however, is not known. Additionally, FISH was 
a useful tool in this study as G. sulfurreducens was 
highly enriched and could be easily targeted. 

The clone library analysis was the only method 
that showed a possible change in the microbial 
community for the underperforming replicate 
UAJA3, indicating a reduced predominance of 
G. sulfurreducens and a greater predominance of 
clones (36%) most similar to Bacteroidetes, but at a 
very low (89%) identity. This change in the micro- 
bial community for the UAJA3 replicate does not 
necessarily indicate that the microbial community 
was the cause of the reduced performance of this 
replicate. It is possible that the performance reflects 
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Table 3 Primer hits to the major phylogenetic taxa found in this study using RDP probe match 



Molecular Primer Sequence (5' to 3') References RDP probe match values (%) 
technique set 

Deltaproteobacteria Geobacter Bacteroidia 



Clone 


a 530F/ 


GTCCCAGCMGCCGCGG 


This study, Wani et al, 2006, 


0.1 


0 


0.1 


libraries 


1494R 


GGTTACCTTGTTACCACTT 


0.3 


0.2 


0 




530F/ 


GTGCCAGCMGCCGCGG 


Madrid et ah, 2001 


97 


97 


97 




1494R 


GGYTACCTTGTTACCACTT 




0.3 


0.2 


0 




27F/ 


AGAGTTTGATCCTGGCTCAG 


Loffler et al., 2000; Arturrson and 


26 


28 


24 




1541R 


AAGGAGGTGATCCAGCC 


Jansson, 2003; Kiely et al, 2011b 


3.8 


6 


0.3 




63/ 


CAGGCCTAACACATGCAAGTC 


Lucas et al, 2003, 


0.8 


0.2 


6 




1389 


ACGGGCGGTGTGTACAAG 


Witt et al., 2011 


78 


82 


31 




41F/ 


GCTCAGATTGAACGCTGGCG 


Hongoh et al., 2003; 


6.4 


0.2 


0 




1389R 


ACGGGCGGTGTGTACAAG 


Fisher et al., 2007 


78 


82 


31 


DGGE 


341F/ 


CCTACGGGAGGCAGCAG 


Watanabe et al., 2001; 


96 


98 


98 




534R 


ATTACCGCGGCTGCTGG 


Xing et al, 2006 


86 


98 


98 




341F/ 


CCTACGGGAGGCAGCAG 


Watanabe et al, 2001; 


96 


98 


98 




926R 


CCGTCAATTCMTTTGAGTTT 


Xing et al, 2006 


80 


96 


96 




968F/ 


AACGCGAAGAACCTTAC 


This study, 


59 


15 


0.5 




1401R 


CGGTGTGTACAAGACCC 


Watanabe et al, 2001; 


1 


1 


0.3 








Xing et al, 2006 








454 


343F/ 


TACGGRAGGCAGCAG 


This study, 


96 


98 


98 


Sequencing 


926R 


CCGTCAATTYYTTTRAGTTT 


Liu et al, 2007 


96 


96 


96 






Wang and Qian, 2009 










27F/ 


AGAGTTTGATCCTGGCTCAG 


Kunin et al, 2010, 


27 


28 


24 




342R 


CTGCTGCSYCCCGTAG 




96 


98 


98 




1114F/ 


GCAACGAGCGCAACCC 


Engelbrektson et al, 2010; 


94 


97 


0.7 




1392R 


ACGGGCGGTGTGTRC 


Kunin et al, 2010 


80 


86 


37 


FISH 


GBS207 


GAAGACAGGAGGCCCGAAA 


Richter et al, 2007 


0.1 


3 


0 



Abbreviations: DGGE, denaturing gradient gel electrophoresis; FISH, fluorescent in-situ hybridization, RDP, ribosomal database project. 
a This primer sequence was incorrectly written in the reference. The correct sequence is S'-GTGCCAGCMGCCGCGG-S' (Madrid et al, 2001; 
Rossetti et al, 2003). Deltaproteobacteria, Geobacter and Bacteroidia all had 97% probe match values when the correct sequence was used. The 
specificity of the PCR conditions would determine whether one mismatch would allow a primer to anneal to the template. This one mismatch did 
not appear to bias the targeting of any particular taxa as all probe match values for the incorrect primer gave ~0% match. 



an undetectable physical difference of the reactor, 
such as a poor seal with a gasket allowing additional 
oxygen to leak into the reactor. Even this possible 
difference in the UAJA3 community cannot be 
certain because of the lack of agreement between 
the pyrosequencing and clone library approaches for 
this sample. 

Conclusions 

On the basis of performance and microbial commu- 
nity analyses of triplicate MFCs inoculated from 
three different sources, it was determined that 
electrochemical performance of the reactors con- 
verged after 2 months of operation (20 cycles). The 
use of at least two replicates was sufficient to 
characterize these systems in terms of electroche- 
mical performance. When replicates produced simi- 
lar voltages and power densities, community 
analysis of replicates produced similar results. 
Results of this study also suggest that at least one 
non-PCR based technique, such as FISH, should be 
used to verify findings of the community analysis 
and corroborate findings on microbial abundances 
in the biofilms. 
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